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We present a search for neutral supersymmetric Higgs bosons decaying to t pairs produced in 
pp collisions at ^/s = 1.96 TeV. The data, corresponding to 310 pb~^ integrated luminosity, were 
collected with the CDF II detector at the Tevatron collider at Fermilab. No significant excess above 
the standard model backgrounds is observed. We set exclusion limits on the production cross section 
times branching fraction to r pairs for Higgs masses in the range from 90 to 250 GeV/c'^. 

PACS numbers: 14.80.Cp, 13.85.Rm, 12.60.Fr, 12.60.Jv 



One of the outstanding questions in particle physics is 
the dynamics of electroweak (EW) symmetry breaking 
and the origin of particle masses. In the standard model 
(SM), EW symmetry is spontaneously broken through 
the Higgs mechanism , which predicts the existence of 
a massive scalar Higgs boson hsM- Theoretical difficul- 
ties related to divergent radiative corrections to the hsM 



mass have natural solutions in supersymmetric (SUSY) 
models 

The minimal supersymmetric extension of the stan- 
dard model (MSSM) % is the simplest reahstic SUSY 
theory. The Higgs sector in the MSSM consists of two 
charged and three neutral scalar bosons. Assuming CP- 
invariance, one of the neutral bosons {A) is CP-odd, and 
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the other two {h,H) are CP-even. Throughout this Let- 
ter, we use h (H) for the hghter (heavier) CP-even neu- 
tral Higgs boson, and to denote any of h,H,A. At tree 
level, the MSSM Higgs bosons are described by the mass 
of A (m^), and tan/3 = Vu/vd, where w^, Vd are the vac- 
uum expectation values of the neutral Higgs fields that 
couple to up-type and down-type fermions, respectively. 
The Yukawa couplings of A to down-type fermions (such 
as the b quark and r) are enhanced by a factor of tan (3 
relative to the SM. For large tan/? one of the CP-even 
bosons is nearly mass-degenerate with A and has simi- 
lar couplings. The dominant production mechanisms of 
neutral MSSM Higgs bosons at hadron colliders are gluon 
fusion 4] and hh fusion |^ 0| . The leading decay modes 
of A and the corresponding mass-degenerate CP-even 
Higgs boson are 4> ^ hi {^^ 90%) and </) ^ rr (- 10%). 

In this Letter we present the results of a search for 
neutral MSSM Higgs bosons produced in pp collisions at 
■y/s=1.96 TeV. The data sample was collected with the 
CDF H detector at the Fermilab Tevatron between 2002 
and 2004 and has an integrated luminosity of 310 pb~^. 
The search is performed in the <j) ^ tt decay channel 
for 90 < mA < 250 GeV/c^. One r is detected in the 
decay to an e or /i and neutrinos, and the other in the 
decay to hadrons and a neutrino. In the following, we 
use Te, r^, and Thad as shorthand notations for the decay 
modes t eve^r, t — > fJ-J^ni^T, and r — > hadrons Vr, 
respectively. Previous and related searches in the di-r 
channel are presented in 0, ■ 

CDF II is a general purpose detector with tracking 
and calorimetry. The tracking system consists of silicon 
micro-strip detectors and a cylindrical wire drift cham- 
ber. It is immersed in a 1.4 T magnetic field produced 
by a superconducting solenoid. Electromagnetic (EM) 
and hadronic (HAD) sampling calorimeters are located 
outside the solenoid and cover detector pseudorapidity 
I77I < 3.6, where f] — — ln(tan(6'/2)) and 8 is the polar an- 
gle with respect to the proton beam. The calorimeters are 
divided into towers with projective geometry. A central 
electromagnetic shower maximum detector (CES) con- 
sisting of proportional chambers with anode wires paral- 
lel to the beam axis and orthogonal cathode strips is em- 
bedded in the EM calorimeter at a depth of six radiation 
lengths. The CES is used to determine the position of 
EM showers with spatial resolution of ^ 0.5 cm. Muons 
are identified by a system of drift chambers located out- 
side the calorimeter volume with combined coverage ex- 
tending to \r]\ < 1.5. The luminosity is measured by gas 
Cherenkov counters located in the detector forward and 
backward regions (3.7 < < 4.7) with 6% precision pH ). 

The search for cj) tt requires detection of an e or 
fj, (from Te, T^) and the reconstruction of the Thad de- 
cay products. Events are pre-selected with "lepton plus 
track" triggers 0|. The triggers require a lepton (e,/i) 
candidate and another track, both pointing to the cen- 
tral calorimeter (|ry| < 1.0) and having azimuthal sepa- 



ration A(/j > 10°. The overall trigger efficiency for signal 
events passing the selection criteria described below is 
greater than 90%. The algorithms for e and ^ identi- 
fication are described in detail in 0. The vector sum 
of the transverse momenta |l2| of the neutrinos from r 
decays appears as missing transverse energy {fr), de- 
termined from the imbalance of energy deposition in the 
calorimeter 0| . The decay products in Thad form narrow 
jets with low multiplicity of neutral and charged parti- 
cles. The positions and energies of 7r°'s and photons are 
reconstructed with the CES detector and EM calorime- 
ter, respectively. In this search we do not distinguish 
reconstructed photons and 7r°'s, and all neutrals are as- 
sumed to be TT^'s. We use a variable-size signal cone and 
an isolation annulus to identify Thad candidates. The 
four-momentum of Thad is calculated from tracks and tt'^'s 
contained in the signal cone. Particles in the isolation an- 
nulus are used to impose requirements that discriminate 
against quark and gluon jets: the scalar sum of the pt of 
tracks (sum of Et of tt'^'s) in the isolation annulus is re- 
quired to be less than 1 GeV/c (1 GcV). We select Thad 
candidates with one or three tracks in the signal cone 
{Nllg =1,3) with Pt > 1 GeV/c , consistent with the 
dominant r decay modes. In the Nll^ — 3 case the sum 
of the electric charges must be equal to ±1. The invariant 
mass of the hadronic system is required to be less than 
1.8 GeV/c^. Electrons are rejected by imposing the con- 
dition (i;'=7P*[j=)(0.95-/) > 0.1, where E"^ is the energy 
in the calorimeter cluster produced by the Thad candidate, 
/ is the ratio of electromagnetic to hadronic energy in 
the cluster, and Pllg is the scalar sum of track momenta 
in the signal cone. Muons are suppressed by requiring 
E^ > 15 GeV. The Thad identification efficiency increases 
from 38% at transverse momentum of the hadronic sys- 
tem plf'^ = 15 GeV/c to ~ 46% for pl^<"^ > 25 GeV/c. 
The probability for misidentifying a quark or gluon jet as 
Thad is measured using jet data samples. It is ~ 1.5% for 
jet transverse energy Ei^*'—20 GeV, dropping to ^ 0.1% 
for £;^''*=100 GeV. 

The acceptances for signal and most of the back- 
grounds are determined from samples of Monte Carlo 
(MC) simulated events produced by the PYTHIA event 
generator with CTEQ5L |0| parton distribution 
functions (PDF's). Tau decays are simulated by the 
TAUOLA package Detector response is simulated 

with a GEANT-based model of the detector. 

The dominant (and irreducible) background in the final 
sample of selected events is from inclusive Z /^* produc- 
tion with subsequent decays to r pairs. It is estimated us- 
ing MC simulated events with normalization correspond- 
ing to cf{pp Z/-f*) X BR{Z/-f* II) = 254.9 pb in 
the di-lepton mass region 66 < mu < 116 GeV/c^ [l8j |. 
The second largest background contribution comes from 
processes with quark or gluon jets misidentified as Thad, 
such as di-jet and multi-jet, W-|-jets and 7-l-jets produc- 
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tion. These backgrounds are estimated from the data 
by applying jet^ Ttad misidentification rates to jets in 
events that pass all selection criteria except for Tfiad 
identification. The validity of the predictions is verified 
using independent data samples representing the back- 
ground processes. The third group of backgrounds in- 
cludes Z/7* [1 = e,n), WW, WZ, ZZ, and tt pro- 
duction. Their contributions are determined from MC 
samples normalized to the theoretical cross sections. 

The events in the T^Thad {T^iThad) channel are selected 
by requiring one e candidate with p^r^^^ > 10 GeV/c, 
and one r^ad candidate with pl^'^ > 15 GeV/c and op- 
posite electric charge. Low-energy multi-jet backgrounds 
are suppressed by rejecting events with | -I- \p^°''^\ + 
\^t\ < 50 GeV. Backgrounds from W + jet events are 
suppressed by imposing a requirement on the relative di- 
rections of the visible r decay products and fr- We 
define a unit vector ( along the bisector of the angle be- 
tween the directions of e (/i) and Thad in the transverse 
plane. The projections p^'" = ( p'e(M) + ~fhad) ' C and 
= ■ C. are required to satisfy p^^ > 0.6p^^^ — 
10 GeV/c. This condition removes ~ 85% of the W + 
jet events passing the other selection criteria, while re- 
taining ~ 95% of the signal. To suppress backgrounds 
from Z II decays with a misidentified lepton, we do 
not accept events with invariant mass of an e (/i) and a 
single-track Thad candidate within 10 GeV/c^ of the Z 
mass. The combined signal acceptance for a Higgs bo- 
son of mass 90 GeV/c^ (250 GeV/c^) in the r^Thad and 
TfiThad channels is 0.8% (2.0%). 

The systematic uncertainties for particle identifica- 
tion efficiency are 3.5% [rhad), 1.3% (e), and 4.6% (^). 
The uncertainties in trigger efficiency for the T^Tfiad and 
Tf^Thad channels are 2.1% and 1.4%, respectively. The 
uncertainty in the determination of backgrounds due to 
jet^ T misidentification is 20%, resulting in 3% effect on 
the total background estimate. The systematic uncer- 
tainty in signal acceptance from event-level cuts is less 
than 2%. The imprecise knowledge of the PDF's intro- 
duces an additional 5.7% uncertainty on signal accep- 
tance [l^ . 

Figure ^ shows the track multiplicity distribution for 
Thad candidates in the data, along with the background 
predictions. The characteristic enhancement in the one- 
and three-track bins clearly shows the contribution from 
events with Thad in the final state. The total number 
of expected events from SM processes after applying all 
selection criteria is Nsm = 496 ± 5(siat) ± 28{sys) ± 
25{lumi). The contributions from Z/j* ~* tt, back- 
grounds with jet—!- T misidentification, and all remaining 
background sources are 405, 75, and 16, respectively. 

We observe 487 events, in agreement with Nsm- To 
probe for possible Higgs signal we perform binned likeli- 
hood fits of the partially reconstructed mass of the di-r 
system {rrims) defined as the invariant mass of the vis- 
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FIG. 1; Track multiplicity for hadronically decaying tau can- 
didates before applying the opposite charge and require- 
ments. 
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FIG. 2: Example fit of the m^is distribution for signal with 
mA=140 GeV/c^. Signal and background normalizations cor- 
respond to the fit results for signal exclusion at 95% CL. 



ible tau decay products and Tf^T- The backgrounds are 
allowed to float within limits set by Gaussian constraints 
corresponding to the systematic uncertainties in trigger 
efficiencies, particle identification, production cross sec- 
tions, PDF's, event cuts, and luminosity measurement. 
Potential differences in rtims shapes between data and the 
MC simulation in different channels are treated as sys- 
tematic uncertainties. We create signal and background 
fTiTjis templates with the MC energy scales shifted from 
the nominal values according to the uncertainties, and 
study the effect on hypothetical cross section measure- 
ments. The deviations from the results obtained with 
the nominal templates are parameterized in terms of 
Higgs mass and input cross section. An example fit for 
rtiA — 140 GeV/c^ is shown in Figure |21 We observe no 
signal evidence for vtia = 90 — 250 GeV/c^, and set exclu- 
sion limits at 95% CL on aijip ^ + X) x BR((^ — + tt) 
as shown in Figure |31 The sensitivity of the limit-setting 
procedure is determined from MC simulations assum- 
ing no signal. The rfims shape uncertainty leads to 
15% (5%) deterioration of the limits for the low (high) 
end of the considered ttla region. The observed limits 
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FIG. 3: Upper limits at 95% CL on Higgs production cross 
section times branching fraction to r pairs. The expected 
limits from the pseudo-experiments are also shown. 



range from 24.4 pb for uia = 90 GeV/c^, to 9.3 pb for 
niA = 140 GeV/c2, to 1.8 pb for ruA = 250 GeV/c^. 

Using the theoretical predictions for the MSSM Higgs 
production and decay to r pairs we interpret the lim- 
its on (t{pP ^ (j) + X) X BR{(f> — > rr) as exclusions of 
parameter regions in the tan/3 vs vriA plane. The cross 
sections are obtained from SM calculations and scaling 
factors (J MSSM Io'sm accounting for the modified Higgs 
couplings [23 . The cross sections for gluon fusion medi- 
ated by a 6-quark loop are calculated with the HIGLU 
program . The corresponding values for bb ^ <j) + X 
are taken from . The scaling factors and BR[d> tt) 
are calculated with the FeynHiggs program [22|- They 
depend on tua, tan/?, the SU{2) gaugino mass parame- 
ter M2, the SUSY mass scale Msusy, the squark mixing 
parameter Xj, the gluino mass rrig, and the Hi ggs mix- 
ing parameter /i. We consider four benchmarks |23l |: the 
^maa; ^^^^ no- mixing scenarios, with /i > and /i < 0. 
The excluded tan f3 vs niA regions are shown in Figure 0] 

The LEP experiments have excluded tua ^ 93 GeV/c^, 
and higher-mass A for small tan/? |2^. Our search is 
complementary, providing sensitivity in the large tan/3 
region. The excluded parameter space in the tan/? vs 
•niA plane for /.t < is similar to the D0 results obtained 
in the 4> bb decay mode and extends to higher 
mA- Moreover, our results in the (j) —f tt channel allow 
us to set comparable exclusions for scenarios with /i > 0, 
as the lower production cross sections are compensated 
by an increase in BR{(f) — > tt). 
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